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A POSSIBLE SHOCK EFFECT ASSOCIATED WITH SEAQUAKES 


INTRODUCTION 

The term 'seaquake' has been defined to mean simply an earthquake whose 
epicenter lies at sea rather than on land. Initially, the word was used to de- 
scribe any seismic disturbance of the ocean including tsunamis, but the term 
has now been restricted to refer to only the effects of earthquakes felt on board 
vessels at sea (Richter, 1958). While there is little difference between the 
cause of an earthquake at sea and that of any comparable continental earthquake, 
there exist considerable differences in effect. 

Descriptions of seaquakes have remarkable similarities. The first such 
were presented by Rudolph in his papers "Ueber Submarine Erdbeben und 
Eruptionen" published in 1887 and 1895, in which he gives over 400 descriptions 
from ships' logs. At the time he performed his studies, comparatively few of 
the shocks felt on board ship had been recorded by seismic instruments on land. 

It was originally thought that these shocks were probably due to small earth- 
quakes. Improved instrumentation has since shown that most of the earthquakes 
were probably of moderate size (Gutenberg and Richter, 1954). 

In the ship reports of seaquakes collected by Rudolph, although many are of 
a minor nature, there are those that are not minor. A seaquake experienced by 
several ships near the Cape Verde Islands in 1893 was described as "a succession 
of heavy shocks felt, cracking sounds, ship trembled most violently. Vibration 
terrible. " Another was described by Captain John Tadot of 'H.M.S. Sylvia' in 
the North Pacific near Japan in 1870 as "Three distinct shocks of an earthquake. 
Three bump3 as if the ship had grounded and had been lifted up and left again by 
the sea." Captain Garden of the 'Northern Monarch' in 1878 at a position in the 
Indian Ocean 12°4' S, 84° 38' E wrote: "Observed sea thrown up to a great 
height possibly 80 feet or more in a column; this occurred three or four times, 
each upheaval lower than the preceding one; effect similar to that produced by 
a torpedo. Examining the position immediately from aloft, the water at the 
place was observed to break three or four times like heavy breakers and then 
smooth down and we saw no more. The place from the ship about S.E. by E. 
distant 5 to 6 miles . " 

Some earthquakes have produced measured results that are almost unbe- 
lievable. In the instance of the Kwanto Earthquake which occurred on Septem- 
ber 1, 1923 in Japan, a hydrographic survey of Sagami Bay shortly afterward 
revealed permanent peak-to-peak vertical changes in the floor of the bay amount- 
ing to 1500 feet, (Heck, 1965). 
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Descriptions of ordinary seaquakes are distinguished by the following 
features: 

a. Most report a loud rasping or grating sound sometimes accompanied 
by violent shying as if the ship had hit a hard object and was skidding 
over it. 

b. Some report at least one (Eiby, 1957) and sometimes two or more sharp 
jolts against the hull of the ship. 

c. Many are accompanied by sounds variously described from ’distant 
thunder’ to ’a tremendous and explosive noise. ’ 

d. The energy in some shocks is sufficient to cause large ships at dis- 
tances of over 100km from the epicenter to list from 5 to 10 degrees . 

e. Many ships are in deep water, i.e. , > 1000m depth when the seaquake 
is experienced. 

f . Usually there is no apparent motion of the surface of the water connected 
with the shock. 

Let us consider how the observed effects of seaquakes may be brought about. 
If we regard much of the earth as an elastic solid, the sudden release of energy 
during an earthquake generates two types of elastic waves: compression waves 
or P waves and shear or S waves; the designations P and S refer to the primary 
and secondary arrivals at the surface. The compression wave, being faster, 
arrives first and is transmitted from the mantle, through the crust to the floor 
of the ocean. Since water is a much less dense medium than rock, the wave is 
refracted upward and would account for the description of a jolt to a ship. The 
S wave does not propagate through the water, but rather is reflected at the rock/ 
ocean interface and is partially converted to a compression wave. This would 
then propagate upward and would account for a second jolt to a ship of less in- 
tensity than the first and distinctty separated in time from it. Additional shocks 
are caused by secondary P waves , reflected waves , etc . Sounds accompanying 
the seaquakes have a variety of causes. Some are the direct result of the shock 
waves striking the ship. Others can be attributed to the transfer of the shock 
into the air; e.g. the sound of distant thunder can be explained by the refraction 
of sound waves in atmospheric temperature gradients. 

With the exception of the studies performed by Rudolph, effects of seaquakes 
have here-to-fore received little attention. This may be due to the lack of sig- 
nificant damage done to shipping or to the fact that ’in situ’ measurements are 
not commonplace. In either case, of the ships that have experienced a large 
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seaquake, most have been well outside (> 100km) the epic antral area. The 
relationship between earthquake belts and shipping density is given in Figure 1. 


DATA 

On April 29, 1970 at 1401 GMT a large eaithquake of shallow focus and 
registering 7. 5 on the Richter scale- occurred in the Guatemala Basin off the 
coast of the Mexican sh e of Chiapas (CFSLP, NOAA). The location of the epi- 
center was 13. 5°N, 92.o°W. The main shock was preceded by a strong fore- 
shock at 1122 GMT, which reached a magnitude of 6.5R. Between 1122 GMT 
and 2015 GMT over 70 seismic events were recorded; there were no reports of 
damage nor of fissuring. The mareographic station at Salina Cruz, Oaxaca, 
was out of order and no report is available on the tide. 

At 1140 GMT, approximately 18 minutes after the initial foreshock, the 
11. 5ju IR radiometer on board the ITOS-1 spacecraft recorded the thermal in- 
formation presented in Figure 2. A circular area approximately 60km in dia- 
meter in the immediate vicinity of the earthquake epicenter shows an anomalous 
temperature enhancement of approximately +3°K over the surrounding environ- 
ment. An average water temperature of 302°K is indicated by the IR data for 
this area. It is of some interest to note that this is slightly lower than that (86° F 
or 303°K) given by charts of mean sea surface temperature (Figure 3) for April 
1970 prepared by the U.S. Navy Oceanographic Office. Other measurements of 
sea surface temperature by the radiometer at this time, however, are somewhat 
lower still than the corresponding 'in situ* measurements. This is due to the 
radiative transfer of the atmosphere, the effect of which can be determined from 
the radiance. 

The radiance (I) as measured by the s/c radiometer can be approximated 
by the equation 


I = r s B(T S ) +(l-r s )BT (Kunde) 

which assumes a unit level atmosphere where 

B = Planck intensity [@ 11. 5ju] 

T s = Surface Temperature 

t s = Transmissivity of atmosphere 

T = Atmospheric Temperature 
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For the case T s = 303' K, T = 290° K, r s =0.6 the temperature as measured 
by the spacecraft sensor would be 298° K or about 5°K lower than the actual sur- 
face temperature. The presence of any clouds in the field of view of the sensor 
would serve to lower the radiance even more. 

The 3°K temperature anomaly is not observable in ITOS-1 or Nimbus 4 data 
in other orbits prior to or following the earthquake. 

On August 11, 1970 at 1020 GMT a major shallow focus earthquake of 
R’chter magnitude 7.6 occurred in the region of the New Herbrides Islands in 
the South Pacific (CFSLP, NOAA). At 1255 GMT the area was observed by the 
11. 5 m IR radiometer on board the Nimbus 4 spacecraft. The water temperature 
in the immediate vicinity of the earthquake, i.e. , 14. 1°S, 166. 7° E, was re- 
corded to be between 296° and 297 °K or about 2° warmer than the surrounding 
area. 

In both the case of the Guatemala Basin seaquake and that of the New 
Hebrides quake it is postulated that the observed temperature effects could be 
the direct result of a shock wave caused by Hie quake and propagating through 

water and atmosphere. y 

DISCUSSION 

The {shock, sound, elastic, c -mpression} wave resulting from the release 
of (stress, strain} energy by an earthquake whose hypocenter lies in the upper 
mantle travels through the upper mant.-: at speeds of from 7 to 8 km/sec 
(Richter, 1958). Between the mantle and the crust is a sharp, almost world- 
wide discontinuity known as the Mohorovicic discontinuity or 'Moho\ Above 

this is the crist, a layered structure of decreasing seismic velocities, the r 

uppermost portion of which is composed of sediments with very low velocities 
and very high attenuation (Stacey, 1969). In instances where it has been meas- 
ured, there is a direct correlation between depth and displacement. In the case 
of the Idu (Japan) earthquake, the measured displacement at the surface was 
three feet, whereas in the tunnel of Toma on the Tokyo-Kobe line 530 feet below, 
the displacement was 8 feet (Tazieff, 1964). 

The average thickness of the continental crust is between 30 and 40km while 
the average thickness of the oceanic crust is only 5km and of much simpler and \ 

more uniform basaltic structure (Stacey, 1969). A shock wave being propagated ; 

upward from a continental hypocenter at a depth of 30km would be severely at- 
tenuated and fractionated by the crustal structure and in particular by the in- 
homogeneous fractures near the surface. In the case of a seaquake, propagation 
is through a thin crust followed by the homogeneous medium of water. Thus, 
less attenuation would be expected. 
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The compression wave from the earthquake propagates as a shock wave 
through the water. Since the piston motion of the ocean floor is finite, the shock 
wave must he followed by a rarefaction. 

The rarefaction wave velocity is the sound velocity behind the shock plus 
the particle velocity (u). A rarefaction wave will therefore overtake the shock 
front and, by cooling and expanding the fluid behind the front, will cause a rapid 
weakening of the shock. It can be readily shown, then, that shock propagation is 
a critical function of the relationship between the piston like motion of the oc'" >.n 
floor and the depth. If the initial compression wave has a duration of 0. 1 sec 
the rarefaction wave would overtake it at a depth of less than 500 meters . (This 
is illustrated by the characteristic diagram for a 6 kilobar shock given in Fig- 
ure 5.) If the total water depth is ~ 2500 meters, as in the case of the Guatemala 
Basin earthquake, any effect upon the surface would be small for reasonable 
displacements of the ocean floor. A value t > 0.25 sec is required to produce 
the effect seen in the IR data. This would require a peak deflection of the ocean 
bottom over the epicenter in excess of 70 meters. The water at the surface 
would not necessarily break since the height/v.'idth (60km) ratio would only be 
on the order of 0. 001; however, the effect upon an object in the immediate vi- 
cinity of the center of such a shock would be catastrophic. * As one moves away 
from the epicenter, the shock strength would rapidly diminish due to both the 
dissipation of the shock front and the effects of rarefaction waves. 

When a shock wave moves with a velocity U through a stationary fluid of 
initial pressure P 0 and specific volume v 0 , the pressure P, specific volume v, 
and particle velocity u, of the fluid behind the shock are determined by the 
Rankine-Hugoniot conditions (Rankine, 1870, Hugoniot, 1887). These express 
the conservation of mass, energy, and momentum of an element of the fluid 
through which the shock front moves and are given by: 

U = [(p - Po) (v 0 - v)] H 

u = V 0 [(p - Po)/(v 0 - v)] l/j 

AH = (1/2) (p - po) (v + v 0 ) 

where AH is the specific enthalpy (defined as £ internal energy per gram+ pv) 
increment of an element of fluid when it passes through the shock front 
(Richardson, 1947). In addition to the above relationships one also needs the 


♦The characteristic time Op, is the time for a plate to deflect to a fraction 1 - 1/e = 0.63 of its final set in 
response to a step pressure and is called the plastic time by Kirkwood. O p for a diaphram 10.0 inches in 
diameter and 0.1 inch thick with a yield stress of 60,000 p.s.i. is 600jusec, (Cole, 1948). 



speed of sound in the fluid c, the Riemann function a , the undissipated enthalpy 
w, and the equation of state of water. 

c = (dp/dp) A $ p = l./v, s = entropy 

a = f [v(p' ,S)/c(p' ,S)]dp' 

J Po 

<*> = | v(p' , S)dp' . 

Tait equation of state: 

v(0, T)-v(p, T)/v(0, T) = (l/n) log [l + p/B(t)] 

where t = (T-273. 16)°C; n and B(t) are empirically determined. 

The solutions of the above equations which give the desired variables u, U, 
c, v, and T as functions of p involve a series of successive approximations and 
have been carefully dealt with by many authors. A comparison of results is 
given in Table 1 for two pressures, 5kb and lOkb. 


Table 1 

Initial Temperature T 0 = 20° C, P 0 = 1 atmosphere 



P(W» 

u (m/sec) 

U (m/sec) 

c (m/sec) 

v (cc/gm) 

T (°C) 

Richardson 

Hi 

251.0 



0.8668 

— 

Sternberg 

5 

248.8 

2013.7 

2264.5 

0.8781 

32.26* 

Rice 

5 

250.0 

2008.0 

2350.0 

0.8773 

36.0 

Kirkwood 

5 

— 

— 

— 

— 

35.7 

Arons** 

5 

240.2 

2038.9 

— 

— 

— 

Richardson 

10 

425.5 

2335. 0 

2755.0 

0.8120 

— 

Sternberg 

10 

425.1 

2356. 6 

2707.9 

0.8211 

50.48* 

Rice 

10 

426.0 

2352.0 



54.0 

Kirkwood 

10 

— 

— 

— 

— 

55.5 

Arons** 

10 

406.9 

2395. 6 

— 

— 

— 


‘calculated for salt water 
**T o = 25°C 
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Using data supplied by the United States Naval Ordnance Laboratory 
(Sternberg, 1972) for explosions in salt water, the pressure/temperature states 
of a fluid immediately behind a shock front can be plotted on a curve, commonly 
referred to as a "Hugoniot" (Figure 4). An element of fluid initially at a state 
(p 0 , T 0 ) reaches a state of (p , , T ,) immediately behind the shock and then re- 
turns to a state (p 0 , T 0 ' ) along the adiabatic (isentropic) curve. 

In Figure 4, P 0 = 1 bar, T 0 = 20°C. To raise the temperature of the water 
at the surface from T 0 = 20°C so that 22°C < T 0 ' < 23 8 C, would require a sur- 
face pressure p, at the ocean/air interface of between 6 and 8 kilobars, or 
~ 100, OOOp. s. i. (14513 p.s.i. = 1 kilobar). This pressure is probably higher 
than is necessary to account for the thermal anomaly observed by the spacecraft 
sensors since it assumes the full temperature rise is in the water and neglects 
the effect of the propagation of the shock into the atmosphere. 

For a 6 kilobar shock c = 2358.43 msec' 1 and substitution of the correspond- 
ing values of 


p 0 = 1 bar = 1.01355nm' 2 
v 0 = 1. 0134 x 10' 3 m 3 kg' 1 
p, = 6 x 10 8 nm' 2 
v, = 0.864006 x 10' 3 m 3 kg-i 
into the Rankine-Hugoniot equations for u and U give 

u = 287, 55m sec -1 
U = 2090.38 msec' 1 

Figure 5 is a characteristic diagram for a 6 kilobar shock showing the re- 
lationships among the particle velocity (u t ), the velocity of the compressive 
shock (U t ), and the velocities of rarefaction waves associated with the shock. 

The relationship between the total energy and the magnitude of an earthquake 
is given by the empirical equation 

log jo E = 12.24 + (1.44) (m) (Blth, 1966) 

where E is measured in ergs and m is the surface wave magnitude. Substitution 
of the magnitude 6. 5 of the Chiapas shock gives an energy of — 6 x 10 2> ergs. 

The energy needed to heat by 2°C the total volume of water (diameter 60km, 
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depth 10km) indicated by the ITOS-1 temperature anomaly is ~5 x 10 19 or 1% 
of the total energy of the earthquake. 

A 6 kilobar shock propagating through the water gives a particle velocity 
(M 3 ) at the air-sea interface of 288 m/sec. The speed of sound (a,) in air at 
20°C is 344 m/sec (Handbook of Chemistry and Physics, 1959). Using shock 
tables (Liepmann and Roshko, 1957), it is possible to plot the particle velocity 
behind the shock in air as a function of shock Mach number (M) [Figure 6]. It 
is seen that a particle velocity of 288 m/sec corresponds to a shock Mach num- 
ber of 1. 65. 

The entropy generated by a shock is related to the change of total pressure 
£ nd is defined as 


AS d 

— = In i-21 (Liepmann and Roshko, 1957) 
R P 02 

v here p 0 is the stagnation pressure and R is a characteristic gas constant. For 
M = 1.65, 


= 0.376, giving ^21 ■= 1.14. 

Poi P 02 


AS 

which to the first order gives — = 0. 14. 


The temperature change between the initial unshocked gas and the adiabati- 
cally expanded gas after the shock will be small. Therefore, it is permissible 
to set 


AS 


T CgdT _ Cp(Tf -Ti) 
Jti T Ti 


where Cp is the specific heat at constant pressure and T, and Ti are the final 
and initial temperatures, respecii/ely. 


Taking 


TiAS 
R Ti 


= 0. 14 gives 


TiAS = Cp (AT) 


Cp(AT) = (0.14) (RTi) 



where Cv is the specific heat at constant volume and y is the (Cp/Cv) ratio. 


AT 

T 


0. 14 



0.04 


If Ti = 300° K, we obtain AT = 12° C. 


Increasing the atmospheric temperature (T) used in Kunde's equation by 12° 
would have the effect of raising the observed IR surface temperature to a value 
almost equal to the actual surface temperature (T s ), i.e. , 302° K. 


Truttse (1971) has suggested that increased density in the upper atmosphere 
is traceable to the warming of the atmosphere by shock waves generated in the 
troposphere during powerful earthquakes and that most spectacular effects are 
observed when the epicenters are located at shallow depths. Our estimates 
agree with Truttse's hypothesis. 


SUMMARY 

The possibility of cohesive shock waves being propagated through the ocean 
during an earthquake at sea was first suggested by descriptions taken from ships' 
logs published in Rudolphs’ papers almost a century ago. Recent reports are 
similar (Hoffmeister, 1971). A search of infrared photographs taken by two 
NASA spacecraft shortly after the occurrences of some large earthquakes at sea 
have indicated the possibility of a thermal effect associated with the quakes. 

The attribution of this effect to shock waves has been explored. A consistent 
explanation for thermal effects of seaquakes in terms of shock wave propagation 
in water and air has been offered and found to be In quantitative agreement with 
observed results. 
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Figure la. Earthquake Belts 
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Figure 4. Hugoniot and Isentropic Curves for 6 and 8 Id lobar Shocks 







among the Particle Velocity (jit), the Velocity of the Compressive 
Shock (Ut), and the Velocities of the Rarefaction Waves associated 
with the Shock. 
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VELOCITY RATIO 



VELOCITY (rn/sec) 



